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We discuss the partitioning of the Hilbert space of a quantum system induced by the interaction 
with another system at tliermal equihbrium, showing that the higher the temperature the more 
effective is tlie formation of Zeno subspaces. We show that our analysis keeps its vahdity even in 
the case of interaction with a bosonic reservoir, provided appropriate limitations of the relevant 
bandwidth. 
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I. INTRODUCTION 

The quantum Zeno effect (QZE), in its original form, is 
the inhibition of the dynamics of a physical system due 
to frequent measurements of its state [ij. In fact, the 
combined action of unitary evolutions for a short time 
and projection operators gives rise to an effective quan- 
tum evolution which becomes closer and closer to a sim- 
ple projection on the initial state of the system, when the 
time interval between two measurements becomes shorter 
and shorter. This effect has been demonstrated in various 
physical systems 0, 0] ■ 

Subsequent studies have raised the problem of analyz- 
ing the inhibition of the dynamics induced by contin- 
uous measurements, meant as dissipative processes Q- 
Indeed, the emission of radiation from a quantum sys- 
tem could be thought of as measurement process, since 
it says to the observer in which state the system was be- 
fore the emission. Apart from philosophical discussions, 
the presence of a strong decay can be responsible for a 
dynamical decoupling, hindering the dynamics induced 
by other couphngs 

Further developments has brought to the idea that 
even a unitary coupling can be responsible for the inhibi- 
tion of the dynamics induced by smaller couplings [9l-[T]| , 
eventually leading to the concept of Zeno subspaces, i.e., 
the formation of suitable invariant subspaces that can- 
cels the effects of weak couplings, as shown by Facchi and 
Pascazio [T^]. The partitioning of the Hilbert space can 
be obtained also by 'bang-bang' (BB) decoupling. The 
theories of QZE and BB have been recently unified [Tsl - 
[isj . It is worth mentioning that A. Peres has shown how 
similar behaviors can be observed even in the frame of 
classical theory [l^ . 

The importance of QZE is related to conceptual as- 
pects and to the foundation of quantum mechanics p7| . 
but it is also witnessed by a variety of applications in the 
fields of quantum information and nanotechnologies (Tsl - 

US. 

About a decade ago, Ruseckas has studied the influ- 
ence of the temperature of the detector on the quantum 



Zeno effect, showing that a higher temperature of the 
measurement apparatus can enhance the inhibition of the 
dynamics at short time, then resulting in a more evident 
Zeno effect by pulsed measurements 2l|. More recently, 
Maniscalco et al have analyzed the crossover from quan- 
tum Zeno effect to Inverse Zeno effect (IZE, meant as 
the enhancement of the dynamics due to repeated mea- 
surements 0]) in the Quantum Brownian motion, also 
bringing to light a certain role of temperature in the oc- 
currence of such effects [l^] . Very recently, Bhaktavatsala 
and Kurizki have analyzed the influence of the statistics 
of the environment (fermions vs bosons) to the QZE-IZE 
crossover [23j . 

In this paper we study the role of temperature in the 
partitioning of the Hilbert space of a physical system, 
i.e. we analyze the formation of Zeno subspaces due to 
the thermal bosons a system is interacting with. Start- 
ing from the study — analytical and numerical — of the 
interaction between a three-state system and a single har- 
monic oscillator at thermal equilibrium, and between the 
same three-state system and a finite number of thermal- 
ized harmonic oscillators, we find a universal feature cor- 
responding to the partitioning of the Hilbert space of 
the three-state system at high temperatures. Our analy- 
sis allows us to single out a critical temperature for the 
observation of such a partitioning. The scaling of this 
threshold temperature easily allows a generalization to 
the case of infinite harmonic oscillators. 

In the next section we show the basic of this effect by 
studying the single harmonic oscillator case. In the third 
section, we extend the discussion to the case of a finite 
but arbitrary number of harmonic oscillators, and in the 
fourth section we give an estimation of a temperature 
high enough to observe the formation of Zeno subspaces. 
Finally, in the last section we provide a detailed discus- 
sion of the results. 



II. SINGLE HARMONIC OSCILLATOR 
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We consider a three-state system coupled with an har- 
monic oscillator tuned close to the transition 2^3. 
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The relevant Hamiltonian is (fi = 1): 

H = ojk\k){k\+m){'2\ + m{M) 

fc=l,2,3 

+ a;ata + .g(a|2) (3| +at |3) (2|) , (1) 

where |fc) with k = 1,2,3 denotes the generic state of 
the three-state system, ujk is the relevant energy, ft gives 
the strength of the coupling between 1 and 2, a; is the 
frequency of the harmonic oscillator (whose annihilation 
and creation operators are a and a^) and g is the cou- 
pling strength between the three-state system and the 
oscillator. 

For the sake of simplicity wc are assuming that the 
coupling constant g is real, but all the results we shall 
show are valid also in the most general case where g is 
complex and the coupling is ga\2) (3| + g*a) |3) (2| (see 
the discussion in the end of Appendix IB)) . 

This Hamiltonian is structured as a set of invariant 
blocks of the form: 

(cji + nui n \ 

n uj2 + nuj_ gsjn + 1 , (2) 

g\Jn+ 1 W3 + (n + l)w / 

corresponding to the triplet |1) |2) |3) jn + 1). 

The relevant evolution operator is Un{t) = 
exp{—iHnt). When the condition g^/rT+T ^ is 
fulfilled, such unitary operator does not significantly 
change the population of the state |1) |n). Quali- 
tatively, wc can say that in such case the coupling 
57(|1) (2| + |2) (1|) is a small perturbation which docs not 
affect much the dynamics of the state |1) \ti) which is an 
eigenstate in the case fi = 0. 

Suppose now that the two subsystems are initially un- 
corrclated and that the harmonic oscillator is at thermal 
equilibrium: 

P = PA® PB , (3) 

with 

PA = |1) (1| , (4) 



PB = Ypn\n) {n\ , Pn = Z ^exp^--^— ^j, (5) 

where Z =[l- cyi^{-uj / [ksT))]-^ . 
The state at time t can be written as: 

p(t)=^p„C/„(t)|n)|l)(l|(n|f/t(i), (6) 

n 

SO that the survival probability of the atomic state |1) is 
given by: 

i'W-E^'-KilHc/nWI")!!)!'- (7) 

n 



We shall rigorously prove that for any e > it is pos- 
sible to find a temperature such that for T > 
it is P{t) > 1 — e at every time. In the Appendix 
IB] we prove that for any < e < 1 it is possible to 
find an index such that for n > it turns out 
|(1| {n\ Un{t) \n) > "v/l — e for every t (consider the 
special case D = 1, with ci = g-y/rT+T). Therefore, the 
survival probability satisfies the following relations: 

Pit) > 5]p„|(l|(n|C/„(0|n>|l)|^ 

n<nE 

+ E pnVl - e > E p„Vl - e 

= exp(-^)yT^, (8) 

where we have used the relation J2n>n ~ a;"^/(l — a;). 

Let us now consider a temperature T > Tg ^ 
— (2cL>ne)/(fcB log(l — e)), for which one has 
exp(— a;ne/(fcsr£)) > \/l — e. Then, for such high 
temperatures one has P{t) > 1 — e for any t. This is 
a clear manifestation of the formation of temperature- 
induced Zeno subspaces. Indeed, the higher the 
temperature, the more the coupling between atomic 
states |1) and |2) is neutralized by the coupling between 
|2) and |3) mediated by the harmonic oscillator. 

It is important to note that when T « the only block 
really involved in the dynamics is the one corresponding 
to n = 0, and the inhibition of the time evolution can 
occur only ii g ^ U,. On the contrary, when the temper- 
ature increases, there are a number of blocks effectively 
involved in the dynamics, and in most of them it happens 
that the coupling constant g\/n + 1 exceeds fi, even if g 
itself is not large. In the T — oo limit all the blocks (an 
infinite number) are involved and in the majority of them 
the condition g^/n + I is fulfilled. 

To support our previous analysis, we show in Fig. [T] 
the evolution of the survival probability of the state |1) 
when the system is interacting with an oscillator resonant 
to the transition 2 — > 3, for different temperatures. It is 
well visible that as the temperature increases the survival 
probability tends toward unity at every time. 



III. MANY HARMONIC OSCILLATORS 

The previous result can be generalized to a set of D 
harmonic oscillators interacting with a three-state sys- 
tem. The relevant Hamiltonian reads: 

H = ^k\k){k\ + n{\l){2\ + \2){l\) 

fe=l,2,3 
D D 

+ + E 9k{ak |2) (3| + 4 |3) (2|) . 

fc=l k=l 

(9) 
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FIG. 1: (Color online) The survival probability of the atomic 
state |1) as a function of time (in units of at dif- 

ferent temperatures: fcsr/a;23 ~ 0.1 (green solid line), 
kBT/LJ23 ~ 1 (red dashed line), kBT/u;23 = 10 (blue dot- 
ted line), ksT /uj23 = 100 (black bold line). Here uJi/Q. = 20, 
L02/O. = 19, a;3 = 0, L02-3 ^ L02 - LJ3, D = 1, g/Q, = 1. The 
harmonic oscillator is tuned to the 2 <=i 3 transition. 



ni...nn — 



with S = '^i^rikUJk- The only non vanishing terms of 
this matrix belong to the second row, to the second 
column (all the couplings involve state \2) \ni, ...no)) 
and to the diagonal (free energies). Similarly to the 
case of one harmonic oscillator, when it happens that 
gk^/nk + 1 3> Vfc, one can assume that the Hamil- 
tonian in (jlOp with i7 = is the 'unperturbed' Hamil- 
tonian, which has |1) \ni,...nD) as eigenstate; then one 
can consider the effects of the perturbation (the coupling 
ri(|l) (2| -I- |2) (1|)), which does not change much the pre- 
vious situation, leaving the state |1) |ni, n^) as an 
eigenstate up to terms of the order fi/ {gkVnk + !)■ This 
implies that for very large ni, 712, ... njj the population 
of the state |1) \ni, ...no) does not significantly change. 

The circumstance that the state |1) no) is close 

to an eigenstate of the Hamiltonian for large excitation 
numbers is rigorously proven in the Appendix |BJ 

Similarly to the case of a single harmonic oscillator, 
for large enough temperature the blocks with large ex- 
citation numbers are more and more involved in the dy- 
namics and the subspace separation becomes more and 
more significant. 

In Fig. [5] we show the time evolution of the survival 



Each set {|1) |ni, 712, n^), |2) |ni, 712, n^), 
|3) |rti -I- l,n2,...,n_D), ...|3) |ni, 712, n_D -|- 1)} indi- 
vidualizes an invariant subspace, corresponding to the 
following invariant block: 



(10) 



I 

probability of the state |1) when the three-state system 
is interacting with a finite number of oscillators having 
frequencies close to the 2 — > 3 transition. In these calcu- 
lations we have taken g = 0.5 and = 4, in order to keep 
constant the quantity g^D, which resembles the zero- 
temperature decay rate in a bath T = q'^(uj) D(ujV where 
D(uj) is the density of modes of frequency w [2J, [2^. In 
order to better afford numerical calculations in the pres- 
ence of a larger Hilbert space, we have considered values 
of the temperature lower than those considered in the 
case D = 1. The Zeno-like effect is still appreciable any- 
way. Indeed, similarly to the case of a single harmonic 
oscillator, the trend of the survival probability as tem- 
perature increases is well visible. 



IV. THRESHOLD TEMPERATURE 

Let us now study in detail the dependence of the ap- 
pearance of Zeno subspaces on the temperature. 

Generally speaking, if the initial state of a quantum 
system has an overlap y^e*"^ with an eigenstate of the 
Hamiltonian, such that x > 1/2, then the survival proba- 



/ Lui + s n 

n UJ2 + S gWni + 1 

giy/ni + 1 ui^ + S + uJi 



\ goy/no + 1 







go^/no + 1 




UJ3 + 5 + LJo / 
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FIG. 2: (Color online) The survival probability of the atomic 
state |1) as a function of time (in units of fi"^), at different 
temperatures: ksTjixi'iz = 0.1 (green solid Une), kBT/u)23 — 1 
(red dashed line), ksT /uj23 = 5 (blue dotted line), kBT/ui2-i ~ 
10 (black bold line). Here uji/Q = 20, UJ2/O = 19, UJ3 = 0, 
iJ2s = 1^2 —1^3, -D = 4, gfc/fi = g/il = 0.5. The four oscillators 
are tuned to the frequencies: cji = CJ23, ^^2 ~ 0.996 cj23, 1^3 ~ 
0.992 CJ23, ^4 ~ 0.987 CJ23- 



bility can never be smaller than (%— (1 —x))^ = (2% — 1)^. 
We prove this statement in the Appendix |X] 

From Appendix IB] we know that in a given block char- 
acterized by rii, 712, "-D with c = gfn-i large enough, 

there exists an eigenstate, say whose overlap 

with the state 1^ 



ni,n2,...,nr) 



) = |1) |ni,n2, ...,nD) is: 



nl 



> 



n-2,...,nD I 



mc 



(11) 



where m = minj(|a;3 — loi + Ljj\) and A/ = maxj(|cj3 — 

CJl + Wjl). 

To obtain the previous formula, according to our anal- 
ysis in Appendix [B] we require either cuj > wi — \fk 
or ujj < cji — V/c, and 771 > anyway. For 



c > Ce 



4M17 /l + (1 - e)!/-* _ 8^/2^^M 

777 Y 1 - (1 - e)l/4 ^^5^ 



, (12) 



the squared modulus of the overlap is larger than x = 
((1 — eY^^ + l)/2, so that the survival probability of the 



state \A„ 



can never be smaller than v'l — e. 



Let us consider the part of the bosonic Hilbert space 
made of those invariant subspaces where it happens that 
the relevant c is smaller than c^, and introduce the rele- 
vant partition function: 



Ze = zC^gjrij < cl) 



(13) 



The smaller the more the atomic state |1) does not 
evolve, and it turns out that such a state is essentially 
separated from the other two atomic states. 

In the following we will show that the contribution of 
to the total partition function becomes negligible in 
the limit T — !■ cx). 

Calculation on the hypercube — The condition c > 
is surely guaranteed if there is at least one boson num- 
ber 77fc > Tig, where is the smallest integer number 
larger than c^/g^-^^ with g^-m = min^ \gj\. Therefore, the 
condition c > Cg is a fortiori satisfied if 77^ > 77^ Vfc. 

The ratio between and the total partition function 

is: 

^<^^^-f[il-^T), (14) 

^tot ^tot , 

where Z{{nj < n^}) is meant as the partition function 
corresponding to nj < Vj: 



Z{{n, < 77J) 




(15) 



while 



D 

Ztot = n(^"^j)' 



(16) 



exp(— (ijj/(fcBT)). Introducing Wj, 



with 

maxj tjj and x = miiij Xj = exp(— a)max/(fcs2^)), one 
finds 



Z, 



< (1 



(17) 



The request /^tot < 1 — >/! — e = «£, translates into 



T> 



ks log(l - aj'^) 



(18) 



In the limit of small e we have that scales as 



8y2f2M777-ie"i/2 and that log(l 



) « -(e/2)V^ 



provided D is not too large (we need <C 1). There- 
fore, under such hypotheses, it is sufficient to have: 



T > Tc(£) 



(19) 



which eventually guarantees that the survival probability 
of the state |1) is always larger than 1 — e. 

Calculation on the hypersphere — Let us now calculate 
the Zi; partition function using the approximation to the 
continuum, which better applies for high temperature (so 
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that the Boltzmann factors do not change much with 
respect to the n^'s): 



D 



-j:^u,n,/{kBT) -Q 



(20) 



where we have introduced the variables yk — gk\/nk, so 
that duk = 2g^'^ykAyk. 

In the large temperature limit the exponential in 
Eq. (j20p approaches unity and it turns out to be: 



Ze 
Ztot 



2^G 



D c: 



.2D 



(21) 



where Go = {2^Dl)~^ is the appropriate geometric fac- 
tor (see the appendix [C|) , and where we have taken the 
high temperature limit for the total partition function: 

The condition Z^/Ztot < c^e then becomes: 



T > 2G 



2 ~ 

1 /D Ce "-^av 
D 



-1/D 



ks g. 



(22) 



where Pav = (Hnff")^''^ and Wav = {Yln^nY^'^- 

By considering the expression of Cg for small e and that 
of Gd after the Stirling approximation, and taking into 
account that for large D one has (•\/27rD )^/'^ — >■ 1 and 



1, we can just ask: 



T > Ts{s) 



Wav 



m2 gl^ kt 



D 



(23) 



which guarantees that the survival probability of the 
state |1) is always larger than 1 — e. In order to en- 
sure the validity of the previous calculations, condition 
([23|) must always be associated to the high temperature 
limit condition, ksT ^ uJmax^ which allows the passage 
to continuum in Eq. (|20p 

Apart from the technical aspects of the development 
of the calculations (on the hypercube and on the hyper- 
sphere), the two threshold temperatures differ because of 
the regime where they apply and because of the scaling 
with the number of oscillators. In particular, Tc{s) ap- 
plies when the number of oscillators is small, and turns 
out to be proportional to such number. We do not need 
such a threshold temperature to be much higher than 
the frequencies of the oscillators. On the contrary, rs(e) 
applies in the limit of high temperature and can be ap- 
plied even in the case of large number of oscillators. Such 
threshold temperature scales as the inverse of the number 
of oscillators. This means that by assuming that the cou- 
pling constants scale as gt D~^/^, one obtains that the 
threshold temperature Ts{e) is more or less independent 
of the number of oscillators. 




(a) 




(b) 



FIG. 3: (Color online). The two possible physical scenarios 
corresponding to a bandwidth which does not intersect the 
Bohr frequency associated to the 1 <=i 3 transitions. In (a) all 
the frequencies of the oscillators exceed the frequency uji—uj^, 
while in (b) they are always smaller. The numerical simula- 
tions in Fig. [1] and Fig. [2] are both obtained in the scenario 



(b). 



It is worth stressing that, according to the analysis 
reported in the appendix [Bl the previous treatment is 
valid if one of these two conditions is satisfied: either 
a)fc>cji— \/k or uJk<i^i—^3 V/c, always requiring 
m > 0. This means that we have to avoid those situations 
wherein some frequencies of the harmonic oscillators are 
smaller than wi — while some other are larger. In other 
words, the Bohr frequency of the transition 1^3 should 
not belong to the frequency band of the oscillators. The 
two possible scenarios are pictured in Fig. [3] 

Finally, we emphasize that usually Tc(e) significantly 
overestimates the temperature necessary to have the for- 
mation of Zeno subspaces. Moreover, both T > Tc{e) 
and T > Ts{e) provide only suflacient conditions, so that 
in principle the Zeno subspaces can appear for lower tem- 
peratures. 



V. DISCUSSION 

In this paper we have analyzed a simple interaction 
model between a set of harmonic oscillators initially at 
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thermal equilibrium and a three-state system. The latter 
system is prepared in a state (1) which is not directly 
involved in boson-assisted transitions, but it is coupled to 
a state (2) which undergoes transitions toward the third 
state (3) due to the interaction with the bosonic part of 
the system. Therefore there is a coupling between states 
2 and 3 mediated by the harmonic oscillators responsible 
for transitions between such two states. The states 1 and 
2 are directly coupled, but when the coupling between 2 
and 3 mediated by the oscillators becomes very strong, 
a partitioning of the Hilbert space neutralizes the direct 
coupling between 1 and 2. 

There are two ways to make the coupling between 2 
and 3 stronger: one possibility is to consider higher cou- 
pling constants (jfc's, which makes the effect occur even 
at zero temperature; another possibility is to increase the 
temperature of the harmonic oscillators in order to make 
that those bosonic states which stronger couple the tran- 
sitions 2 3 (due to an high value of gk^/nu + 1) are 
much involved in the dynamics. The latter case is the 
one considered in this paper. 

A very remarkable result is the scaling of the threshold 
temperature at which the Zcno subspaces surely appear. 
Indeed, in the limit of high temperature such threshold 
temperature scales as g~^D^^. Therefore, assuming a 
larger and larger number of oscillators, keeping constant 
the quantity g'^^D, we can rigorously assert that even in 
the presence of a bosonic bath the partitioning of the 
Hilbert space occurs. Instead, if the coupling constants 
5fc's scale in such a way that g^.^D significantly decreases 
when D increases, then the temperature necessary to cre- 
ate the Zeno subspaces becomes higher and higher and 
could become infinite when the number of modes of the 
bosonic field is infinite. 

Concerning the limitations of our analysis, it is impor- 
tant to clarify that the proof given in the appendix [B] is 
based on a precise hypothesis about the diagonal entries 
of the matrix. Such hypothesis translates into the as- 
sumption that the frequencies of the oscillators are local- 
ized in a finite band, so that it is possible to find a finite 
maximum and a non-vanishing minimum — both of the 
same sign — for the Bohr frequencies related to transi- 
tions from state 1 to state 3 assisted by the acquisition of 
a bosonic quantum. In other words, the Bohr frequency 
LOi — ws should not be included in the frequency band 
that contains all the frequencies of the harmonic oscilla- 
tors. Therefore, the D ^ oo limit should be considered 
keeping finite the bandwidth of the reservoir and ensur- 
ing that such band does not intersect the Bohr frequency 

Wi - UJ3. 

As a final remark, we mention that if the subsystem 
B is made of spins instead of harmonic oscillators, the 
previous effect seemingly does not occur. Indeed, for the 
partitioning to happen we need to have gk^Jnk + 1 ^ 
which is fulfilled for large n^, and then in the majority 
of the invariant subspaces. Since such large numbers of 
excitations are not possible for a single spin, we guess 
that it is more difficult to obtain Zeno subspaces through 
the interaction with spins at thermal equilibrium, at least 
assuming a structure of interaction analogous to the one 



we considered in this paper. 



Appendix A 

Let be the initial state of the physical system, 

\4>k) the eigenstates of the Hamiltonian. Introducing 
('0(O)|0fc) = ^/Xfce*"*, one can write: 



k 

which evolves into the state 

l'/'(0)> = ^ VHe- 



(Al) 



(A2) 



If there is a Xfe, say Xii which is larger than the sum of 
all the remaining ones, then it is easy to convince oneself 
of the following relations: 



m)m))\ = 1x1 + 

Ixi -^Xk\^ |2xi 



> 



,-i[afc(t)-ai(t)]| 

1|, (A3) 



where we have used J2k Xfc = 1 and we have considered 
the 'worst case' corresponding to ak{t) — ai(t) = tt 
Vfc ^ 1. 

The condition xi > Sfe/i Xfc is clearly equivalent to 
Xi > 1/2. The survival probability is the squared mod- 
ulus of the overlap of the initial state and the evolved 
state: 



(^(0)|V'(i))P> (2x1-1) 



Appendix B 



(A4) 



Let us consider an Hamiltonian of the following form: 



H = 





n 


. 


. 


n 


62 


Cl . 


CD 





Cl 




. 


V 


CD 


. 


■ Sd+2 



(Bl) 



in the basis |^), \B), |Ci), ...\Cd)- Let us assume Sk ^ 
Vfc and introduce the quantities c = (X^feLi c|)^^^, m = 
minf^i(|4+2 - Si\) and M = max£^i(|4+2 - Si\). For 
the sake of simplicity, we are considering the Cfc's as real 
quantities. 

In this appendix we prove that in the limit of very 
large c there is an eigenstate which approaches the state 
1^), provided m > 0. In particular, we will se that there 
exists an eigenvalue Ai such that |Ai — Si\ < 2MU^/c^, 
and that the corresponding eigenstate |Ai) has an overlap 
with \A) which becomes closer and closer to unity as c 
increases. 
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The secular equation — The determinant of the matrix above can be developed (with respect to the first row) as 
follows: 



det H = Si 



( ^2 Ci 

ci h 



CO \ 





- n 



/ n ci 

(53 



\ C_D ... (5d+2 / 



CD \ 





V ... dD+2 J 



where the second minor can be easily evaluated as Jlfet^j^ '^fc- while the first minor is calculated as follows: 



. 




( 


. 


. 




S5 ■ 


. 


- Cl 


. 


■ Sd+2 / 


V 



C2 5i 
C3 5^ 

V ... 5d+2 J V 

£1+2 D+2 

= x{h-c\x{5u~ci>, n + 

k=2 fc=4 k=3,,bfi...D+2 

Taking into account these results, one obtains: 

/D+2 \ / 



:: I \ 




/ Cl 


S3 


. 








C2 





. 







+ C2 


C3 





S5 ■ 


. 


•■ 5d+2 / 




\ CD 







• Sd+2 / 



/D+2 



a-=2 



2^2 



1=1 



S2SI+2 



1=1 



S2S1+2 



The eigenvalue equation, det(_ff — AI) = 0, is obtained through the replacement Sk 



P(A) 

Considering X = 61 + a, one gets 
P{Si+a) = -a 



D+2 

fc=i 



(Si - A) iS2 - A) ^ iS2 - A) iSi+2 - A) 



'D+2 



n -Si-a) 



S2 - Si - a ^ 

rr ^ — ^ 



.fc=3 



1=1 



Si+2 - Si - a 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



The eigenvalue close to 61 
case where Sk+2 > Si Vfc > 1. 
In the limit a -+ one obtains 



Let us first consider the always have to satisfy the previously considered condition 



P{Si) = -1^2 



D+2 



fe=3 



(B7) 



from which one immediately finds P{Si) < 0. 

On the other hand, when c is large enough one can 
prove that P{6i + ^ > for ^ = 2Mn'^/c^. Indeed, 
assume > 2MQ'^/m so that S1+2 — Si — ^ > Vfc > 
1. Therefore, since M > S1+2 ~ <5i — C f^nd hence 
Y^i cf/{Si+2 - (5i - > E; cf/M = cVAf, one obtains: 



S2-Si-^- 



D 



fr{ Si+2 - ^1 - ^ 



< S2-S1- 



c 

2M 



(B8) 



Now, ii §2 — Si < Q then we have a negative quantity for 
any c, and therefore +^) > 0. If instead 62 — Si > 0, 
we need to take c > {2(62 — <5i)A/)^/^. Of course we 



> 2Afri^/m to guarantee 5k+2 — <5i — ^ > 0. 
Summarizing, for c > 



ax {l^VSWw, [2 ((52 - 
(5i)A/]i/2} and ^ = 2Am^/c^ we are sure that P((5i+0 > 
0, and we can then assert that for some rj: < rj < ^ = 
2MVi^ /c? there is an eigenvalue Ai = 5i + 1]. 

Eigenstate and overlap — If ^ and Sk+2 ^Si—rj^ 
Vfc, then one can easily find that the corresponding 
eigenstate is: 



\Si+v) = 



E 



ClT] 



^ iSl+2 -Sl-Tj) 



\Ci) 



(B9) 



which approaches \A) in the limit of very large c. Indeed, 
consider first of all that in such limit we have 77 — )■ 0, 
since < < ^ = 2M^ jc? — )• 0, so that the condition 
Si+2 — (5i — 7^ is easily satisfied. Secondly, one finds 
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|q7?| < \2Mn'^ci/2c^\ < 2Mrj2/c ^ 0, and on this basis 
it is immediate to prove that | + f?|^) | — >■ 1. 

Now, one cannot be sure that (5/+2 ^ Si — ^ > m, even 
for small ^, but assuming > AMVt^ /m one has {5i+2 — 
6i — < 2/m and a fortiori ((5;+2 — (5i — 'rj)~^ < 2/m. 
Then, reminding that rj < 2Mrt^/c, one finds: 



{5i+v\A) 



> 



2 D 

1=1 

— 

AAPn^ 



{Si+2 - Si- r)Y 



D 

E 

1=1 

16M2r22 



mc 



(BIO) 



which clearly approaches unity in the limit c 
this reasoning is invalidated if m = 0. 
The condition 



oo. All 



Now, either 82 ~ 5i > 0, so that the quantity on the 
right-hand side above is positive whatever the number c, 
or 82 — 5i < 0, and that quantity can be made positive 
provided c > {2\S2 — i5i|M)^/2^ This immediately implies 
that P{Si) and P{Si+£,') have opposite signs, provided c 
large enough. Therefore, this time we have an eigenvalue 
Ai = ,5i + 77' with -2Af f72/c2 < 77' < 0. The correspond- 



ing eigenstate \6i - 
we derived for \6i 



rj') is close to \A), since the inequality 
- 77) in (jBlOp is vahd also for \di + rj'). 



Here we are not considering those cases in which 81+2 > 
5i holds for some Ts and does not hold for others. 

Finally, we remark that if the c^'s are not real — corre- 
sponding to the interaction X]fc(fffeafe |2) {^+gla\ |3) (2|), 
which generalizes that in Eq. ([9|) — the previous results 
keep holding, provided the replacement cf. — > |c/j|2 every- 
where. 



Appendix C 



(Ell) 



is satisfied when: 



c > 



4^^M^/X /l + y/l + (1 - x)?t7V(16xM"^ 



m V 2(1 -x) 

(B12) 

Assuming m/VL not diverging and x ~^ 1; we can ne- 
glect the term multiplying 1 — x, then requesting just the 
following: 



> 



AMVL 



1-X 



(B13) 



Of course, we also need to have c > 
max{4Mrj2/TO^ [2(^2 - (5i)M]i/2} ~ which includes the 
previously considered > 2MVP' /m — , but for x close 
enough to unity such condition is surely included in the 
condition (|B13[) . 

Generalization and limitations — All the previous dis- 
cussion keeps its validity if 5k+2 < <5i Vfc > 1, which 
means that n^L3'^('^i ~ '^1) appearing in P{5i) — see 
Eq. dHZl) — has the sign (-1)^. Now, define f = 
—2MVi^/c^, and assume > 2MQ,'^/m in order to have 

Sk+2 — Si—C < Vfc, so that Yli^'^i^i — ^i — O maintains 
its previous sign. Since: 



(B14) 

in order to obtain a change of sign for the quantity P{di + 
^) it is sufficient to have: 



(B15) 



The geometric factor Gd in Eq. (pij) comes from the 
right hand-side integral of Eq. pO|) . In particular, ex- 
ploiting to hyperspherical coordinates, such integral re- 
quires integration over the D — 1 spherical angles. In 
general there are D — 2 angles which can run over [0, tt] 
and one that can run over [0, 27r]. In our particular 
case, the spanned intervals are smaller, since the vari- 
ables 7/j = ffjTij can assume only positive values, and 
hence all the angles span [0,7r/2]. 

The spherical coordinates are defined by: 

yi = r cos 01, 

7/2 = r sin(/)i COS02, 

Us = r sin sin (j)2 cos ^3 , 



yD-i = r sm0i sm02--- sm(/)u_2 cos0D_i 

yo = r sin(/)i sin(/i2--- sin(/)D_2 sin(/)D_i , (CI) 

which correspond to the volume element: 



m-2 



D-l 



Ildy^- = n sin^"'-Vfe X dr n d</.. 



. fe=i 



k=l 



(C2) 



that can be straightforwardly derived from the Jacobian 
of the coordinate transformation \dyk/ dxj\, with Xj = 
for j < D — 1 and xd = r. 
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Therefore one obtains: On this basis we find that the geometric factor raising 

from the integration over the angles is: 



Co = JlJsT ^ (C4) 



k=i Exploiting the Stirling approximation [26|, which is 

2D r- r- valid for large D, one eventually gets: 



'2D 



r - r n sin^f^-^'-Vfecosc/^^d^/), 



^0 

1 c?^ 



X 



2D 2D-^{D-1)\ 2^01 



(C3) 
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